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Abstract 

Long-distance bidirectional transport of organelles depends on the motor proteins kinesin and 
dynein. Using quantitative data obtained from a fungal model system, we previously developed 
ASEP-models of bidirectional motion of motors along unipolar microtubules (MTs) near the cell 
ends of the elongated hyphal cells (herein referred as "unipolar section"). However, recent quantita- 
tive live cell imaging in this system has demonstrated that long-range motility of motors and their 
endosomal cargo mainly occurs along extended antipolar microtubule bundles within the central 
part of the cell (herein referred to as "bipolar section"). Dynein and kinesin-3 motors coordinate 
their activity to move early endosomes (EEs) in a bidirectional fashion, with dynein mediating 
retrograde motility along the unipolar section near the cell poles, whereas kinesin-3 is responsible 
for bidirectional motions along the antipolar section. Here we extend our modelling approach to 
simulate bidirectional motility along an antipolar microtubule bundle. In our model, cargos (par- 
ticles) change direction on each MT with a turning rate SI and the MTs are linked to each other at 
the minus ends where particles can hop between MTs with a rate q\ (obstacle-induced switching 
rate) or 02 (end-induced switching rate). By numerical simulations and mean-field approximations, 
we investigate the distribution of particles along the MTs for different overall densities O. We find 
that even if O is low, the system can exhibit shocks in the density profiles near plus and minus 
ends caused by queueing of particles. We also discuss how the switching rates influence the 
type of motor that dominates the active transport in the bundle. 

PACS numbers: 87.10.Mn, 87.10.Hk, 87.16.Wd 
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I. INTRODUCTION 



Spatial redistribution of organelles is of central importance to all eukaryotic cells. Long- 
distance transport involves the activity of molecular motors that move along polymers of 
tubulin dimers, the so-called microtubules (MTs), powered by the hydrolysis of ATP 
Bidirectional transport along MTs involves the opposing motor molecules kinesin and dynein 
where kinesin takes cargos to the plus ends of MTs and dynein takes the cargo to the minus- 
ends Q| . 

Numerous theoretical studies attempt to describe motility behaviour of cargos along a 
single MT or a network of MTs, e.g., see js-5]. One of the simplest and best studied models 



is the asymmetric simple exclusion process (ASEP) defined on a single track; see review 
and references therein. Extended ASEP models for bidirectional transport have made var- 
ious assumptions to avoid collisions between opposite-directed particles. One is assuming 
that particles are binding and unbinding to/from tracks {l-ioj]. Another is assuming that 
the exclusion principle only applies to particles moving in the same direction and presence of 
a motor in the opposite direction modifies the rate at which motors enter into the site 11]. 

J L 

Alternatively, a high direction-change rate can avoid clusters due to collision [12]. Moreover, 
Evans et al introduced another possibility |13j] to avoid collisions by allowing particle inter- 
changes when they meet. Juhasz introduced a two-lane ASEP [lj] with opposite-directed 
particles moving in separate lanes, thereby avoiding collisions between opposite-directed par- 
ticles. The two-lane ASEP is equivalent to a two-species ASEP in some sense as discussed 
in lsj]. We introduced a multilane model in 16] where particles can change protofilaments 
to avoid collisions, taking into account that a single MT consists of 13 protofilaments 
each of which provides a potential track for motors. 

Based on recent advances in live cell imaging techniques and the use of the fungal 
model system Ustilago maydis, in vivo observation of dynein indicates that collision be- 



tweee opposite-drrected motors rarely occur, W |u|. This a„ows adapting the two-lane 
model in [14j to investigate the bidirectional transport of dynein motors on unipolar MTs. 
The adapted two lane model together with a more sophisticated 13-lane model provides an 
explanation for the formation of dynein accumulation at MT plus-ends Q, 20]. This 
dynein accumulation at MT p lus ends is suggested to prevent the cargo - early endosomes 
(EEs) - falling off the MT [l8(. More recent work on U. maydis has shown that the majority 
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of the fungal cell contains antipolar MT bundles and that unipolar MTs are restricted to 

n 

;he cell poles [21[. Early endosomes - the main cargo of dynein and kinesin-3 in hyphal cells 



22 



23J, undergo long-distance bidirectional motility [24]. Interestingly, bidirectional long- 



distance motility of EEs along the bipolar MTs is mainly mediated (dominated) by kinesin-3, 



whereas dynein is mediating retrograde motility of EEs along the unipolar MT 



2l|. That 



paper also shows that EEs travel over the entire length of the MT array and concludes that 



dynein and kinesin-3 cooperate and that cargo can hop between MTs within the bundle [21 ] . 
They also show that during the EE transport, dynein can detach from the cargo and the EE 
continues motility (carried by kinesin-3) after a short pause. Moreover, short pauses of EEs 
before they continue the directed motility are observed at minus ends of MTs by visualizing 
both minus ends of MTs and EEs J2I]. These suggest that EEs may change MTs at MT 
minus ends by altering their active motor types. 



Hopping between tracks has been modelled for unidirectional traffic |25H29| where a single 
track is followed by parallel tracks. The junction between single and parallel tracks allows 
particles on the single track to step into either of the parallel tracks. Here we provide 
a new model that is based on our previous ASEP model but which includes this junction 
mechanism to describe bidirectional motility of cargo along an antipolar microtubule bundle. 
In our model, the lattice is composed of two MTs that are coupled at microtubule minus 
ends and the arrangement of antipolar MT bundle gives unipolar and bipolar sections within 
the bundle. A detailed description of the model is given in Section [TXJ In Section HTT1 we show 
that the distribution of particles along the lattice can exhibit a variety of phases depending 
on the parameters; in Subsection IIII Al the phases in the unipolar section are analyzed while 
in the following subsections the phases in the entire system are considered. In particular, we 
show that even at a low overall density of particles, particles can accumulate at minus ends 
as well as at plus ends. Moreover, we find a novel type of phase where density profiles of one 
type of particles can smoothly connect between low and high density. The density profiles 
for these different phases are well approximated by our mean field analysis. In Section IIV1 
we consider how the MT switching rates q± t 2 affect the contribution of each type of motors 
to the transport of cargos. Finally, in Section |V] we discuss the biological relevance of this 
model. 



3 



II. A LATTICE MODEL WITH ANTIPOLAR BUNDLING 



In this section, we introduce a simple discrete-lattice model with antipolar bundling of 
two MTs (we refer to each MT as a track). The bundle is presented by a lattice of length 
N and contains two tracks of length N — N\ and N2 respectively that overlap in a common 
section of length iV 2 — N\. The two plus ends of MTs are located at two ends of the lattice 
and the other ends of two MTs - the minus ends (marking the end of the overlap) - are in 
the middle; see Figure [TJ The middle section is referred to as the bipolar section. The left 
sections in the bundle are referred to as the unipolar sections which are of a relative length 
x\ and 1 — X2 respectively where X12 '■= Ni^/N. Considering the symmetric organization 



21] . we assume a symmetric lattice with x\ + X2 = 1 in the 



of the MT bundle shown in 
following. 

In the bundle, each track can support bidirectional transport of particles (particularly EEs 
here) driven by opposite-directed motors (kinesin-3 and dynein); these particles are of two 
types: plus- and minus-type. Plus-type particles are driven by plus-directed motors moving 



towards the plus end of the track whereas vice versa for minus-type particles. Reversa 



transport direction along each track can be realized by a brief "tug-of-war" event 



30-33 



of 



between counteracting motors on the particle. As rare collisions between opposite directed 



motion of EEs in U. maydis are observed 



2 if ] which is similar to the dynein transport, 



we assume transport on each track in a particular direction is on a separate lane as in 
the two-lane model 0] for convenience. Thus, the lattice contains four lanes, each of 
which supports one direction and a single type of particles as illustrated in Figure Q] (b). 



As EEs rarely fall off the MTs 18|, |2l| in the U. maydis, we assume the transport system 



is closed, i.e., there is no injection or exit of particles into/out of the lattice at two ends 
of each track. The closeness of the system leads a particle number conservation during the 
transport. Particles once reach the end of tracks will wait until they change direction to 
continue the motility. 

Further assumptions on the model are made on track switching of particles based on 
experiment observation of EE motility around minus ends in 21] . An antipolar bundle of 
tracks enables particles to switch between tracks together with a change of type and without 
reversing direction. This form of track switching has been observed in U. maydis [2l|. Track 
switching may be possible inside the bipolar section as well as at the junction between 
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bipolar and unipolar sections. For simplification, in this paper we assume track switching 
occurs when passing the junctions associated with minus ends of MTs. Due to the bipolar 
organization in the bundle, plus-type particles can only be on the same track when crossing 
the junctions. When a minus-type particle moves from a unipolar to a bipolar section, the 
minus end in the junction may enhance plus-directed motors on the particle to hop onto 
the other track, resulting in a change of the track on which the particle moves and the type 
of the particle. Thus the minus ends can be viewed as "obstacles" for minus-type particles 
when attempting to step forward on the same track. This assumption on rack switching is 
sufficient to allow particles to travel across the entire lattice from one plus end to the other. 

In our model, we identify each location in the lattice by a pair (/, i) together with a 
transport direction (corresponding to a lane in each track I) where I G {1,2} and i G 
{0, ■ • • , iV — 1} denotes the site along the lattice. Particles hop from one location to another 
with a certain transition rate where the possible transitions we consider are listed below. 

Forward. Plus and minus-type particles on each location can move forward to the nearest 
site along the same lane. We assume equal forward rates for plus- and minus-type 
particles (corresponding to same velocities which is suggested from in vivo trans- 
port 19 1) and homogeneous forward rates on each lane except at the junctions 
(1, N 2 ) and (2, Ni — 1) for minus-type particles. For convenience, we let the equal for- 
ward rate p = 1 and other rates in units of p. For the forward rate between junctions 
for minus-type particles, we set the rate to be 1 — q\ where q\ is the obstacle-induced 
switching rate given below. 

Track switching. We allow minus-type particles at junctions (1, N 2 ) and (l,Ni) to step 
forward onto the second track together with a type change with rate q\ G [0, 1] (called 
obstacle-induced switching rate) and q 2 G [0, 1] (called end-induced switching rate) 
respectively. Similarly, minus-type particles on the second track can switch onto the 
first track at the junctions with the same obstacle- induced/end- induced rate. 

Direction change. Plus- and minus- type particles can change directions on the same 
track and the same site (by changes in type/lane) along the entire lattice. We assume 
homogeneous and equal direction-change rates for plus- and minus-type particles which 
are denoted by w. In addition, we define Q := uN which characterizes the overall 
direction-change rate. 



(a) 




FIG. 1: (a) shows a schematic hypha where two MTs with plus ends at cell poles (tip and septum) 
form an antipolar bundle, (b) shows a schematic diagram of bidirectional transport on the bundle 
which is presented as a discrete lattice of length N. Two tracks in the lattice represent two MTs 
and are of length N — Ni and N2 respectively with plus ends at two ends of the lattice. Plus- and 
minus-type particles on each track move on separate lanes with indicated orientation. The forward 
rate on the same track is set to be 1 except at the junctions between the bipolar section and the 
unipolar sections for minus-type particles. Particles can also change directions on the same track 
with rate uj. We assume that track switching only occurs at the junctions between sections for 
minus-type particles and occurs with rate qi^. If qi t % = then the two MTs are uncoupled. 

The forward and obstacle- induced switching rates for minus-type particles at junctions 
(1, N 2 ) and (2, N\ — 1) are chosen to preserve the overall forward rate p — 1. All transitions 
described above are subject to a simple exclusion principle, i.e., there can be at most one 
particle at each location of the lattice, and presence of a particle at one location prevents 
other particles from moving into that location. 

In the time evolution of this stochastic bidirectional transport, the occupancy at each 
location of the lattice for plus- and minus-type particles at time t, T l ±i (t) changes according 
to the above transition rates (assumed to take place independently and instantaneously) from 
an initial state t 1 ± j(0). The exclusion principle ensures that r l ± G {0, 1}. A special case of 
the model is q\ = qi = where the transport on each of the two tracks is independent of the 
other. In this case the process is no longer ergodic. For the ergodic cases, the statistically 
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stationary state will be independent of initial conditions and a symmetric lattice structure 
gives a symmetric distribution of particles on two tracks where the total numbers of particles 
are equal. In order to compare transport properties with different switching rates q% y %, we 
assume an initial condition satisfying that each track possesses an equal number of particles 
and assume that the system has reached a statistically steady state. 

In statistically stationary state, densities of plus- and minus-type particles on each loca- 
tion (I, i) are defined as the mean occupancy of the particles: 

A ■= °l ■= «-) 

where the brackets (•) denotes the ensemble average. These densities are related to the 
overall density 6 of particles expressed as [3^] 



2(N-N 1 + N 2 ) 2(N-N 1 + N 2 ) 

where 0± are the overall densities of plus- and minus-type particles respectively. The overall 
density O is conserved under the time evolution whereas 0± may not be conserved. Mean- 
while, the mean currents (the mean rate of stepping forward within unit time) away from 
the junctions are given by 

4, = fe(l-ri 4±1 )) 

in units of p on the first track. At the junctions between sections, the mean currents for 
minus-type particles are 

J\ N2 = (l - gO«* a (l - <;v 2 -i)> + ?i«at 2 (1 - <iv 2 -i)> 

and 

J-,m = Q2(t\ Ni (1 - T+ iNl -i)) ■ 

Similarly, we can write the expressions for the mean currents on the second track. The net 
current in the unipolar section is zero and at the junction between sections it is balanced: 

J -,N 2 ~ J +,N 2 -1 + J -,N 2 -l- 

The exact analytical solutions of the density profiles in our model with general parameters 
can be very difficult or impossible to find using methods such as matrix production joj or 



Bethe ansatz 



34j and therefore we use mean-field approximations and numerica 



For the continuous-time discrete-state model, we use a Gillespie algorithm 



simulations. 



351 ] to simulate 



the time evolution. The parameters that govern this bidirectional transport in U. maydis 
are not yet fully known from experiments, but we do set some parameters informed by 
known properties of early endosome transport in U. maydis. We use x\ = 1 — x 2 = 0.2, 
as microtubule minus ends in U. maydis hyphae are shown to be approximately uniformly 
distributed in the middle at about 10% in hyphal length away from cell poles where the MT 
minus ends are "almost" absent 21]. For the parameter £7, in vivo experiments in jisl. \\\ 



suggests a range of run length M = 10 ~ 70 fim in a hyphal length of L = 100 fim which 



gives 
stated. 



371 ] Q G [1, 10]. A lattice length of N = 500 is used for simulations unless otherwise 



A. Symmetry of the system 

The standard unidirectional ASEP on a single lane and the two-lane ASEP developed by 
Juhasz [l^l share a common feature - particle-hole symmetry. This feature does not hold in 
the model we introduce here. However, by adapting the forward stepping rate of plus-type 
particles when crossing the junctions between sections (i.e., assuming there is inhomogeneity 
in stepping rates for plus-type particles at locations (1, N 2 — 1) and (2, A^)), the model can 
still exhibit an exact particle-hole symmetry as explained below. Note that particles are of 
two types, either moving towards or away from the plus ends on each track, we divide holes in 
the lattice into two types. Minus-type holes refer to the holes in lanes for plus-type particles 
(i.e., the first and last lanes in the lattice shown in Figured]), as these holes move towards 
the minus end when the corresponding plus-type particles step forward; plus-type holes 
refer to the holes in lanes for minus-type particles. An (obstacle-induced) track switching 
of minus-type particle from the unipolar to bipolar section indicates an (end-induced) track 
switching of a plus-type hole from the bipolar to unipolar section. Similarly, a plus-type 
particle stepping from the bipolar to unipolar section on the same track indicates that a 
minus-type hole steps from the unipolar to bipolar section. Therefore, if we let the forward 
rate for plus-type particles when crossing the junction at locations (1, N 2 — 1), (2, Ni) be p 
and choose p = 1 — q 2 , then the system possess the particle-hole symmetry, i.e., the system 
is identical under the exchange: 

p+*l-q u qi^q 2 , t^^I-t^. 
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The uncoupled system with homogeneous rates for plus-type particles (i.e., q x ,2 = and 
p — 1) is a special case with particle- hole symmetry. For small q±, the system with homoge- 
neous rates for plus-type particles (i.e., p — 1) still have some approximately "symmetric" 
behaviour. Thus we focus on low and intermediate overall densities; information for high 
overall densities can then be partly deduced from this "symmetry" . 



III. STEADY STATE DISTRIBUTION OF PARTICLES AND MEAN-FIELD 
ANALYSIS 



In the lattice of length N, we take a rescaled position variable defined as x = i/N G [0, 1]. 
In the continuum limit where N — > oo, we re-express J l ±{x) as the mean unidirectional 
current to the plus/minus end on each track (I = 1,2) and J l (x) := J\(x) — J_(x) as the 
net (mean) current. Moreover, we re-express the densities on each track as p\x) and a\x) 
for plus- and minus-type particles respectively. The overall density of plus- and minus-type 
particles can then be re-expressed as 

_ jj(p 1 (x) +P 2 (x))dx _ J^(a 1 (x) + a 2 (x))dx 

+ 2{l-x 1 + x 2 ) ' ~ 2{l-x 1 +x 2 ) 
In the spatial symmetric lattice with X\ + X2 = 1, we have p x {x) = p 2 (l — x) and <r 1 (x) = 

<j 2 (l — x) for these well defined regions. Hence we only need to consider the distribution 
of particles on one track, say the first track, and can ignore the track index superscripts 
without ambiguity. 

As the lattice is composed of unipolar and bipolar sections, we consider these sections 
separately with appropriate boundary rates a± b \ ft±) see Figure El The stochastic process 
in each section is an ASEP similar to that discussed in ujl 20 1 but with different boundary 

n n 

conditions. In both ASEPs in 14!, I20J, the mean-field approximation predicts the density 
profiles well for statistically stationary states. For the mean-field approximation, we ig- 
nore two-point correlations. The density profile in the steady state for plus- and minus-type 
particles in the unipolar/bipolar section of the lattice is then governed by the following equa- 
tions (ignoring the second derivative and second order of 1/N) with appropriate boundary 
conditions (see Ij], Q for details) 

= (2p-l)^-Q(p-a) (1) 

. da „ . . 
= (l-2a)— + tt(p-cr). 
dx 



Meanwhile, the unidirectional currents within each section are given by equations 



J+(x) = p(x)(l — p(x)), J-(x) = cr(x)(l — cr(x)) 
and the net current reads as 

J(x) = p(x)(l — p{x)) — cr(x)(l — cr(x)). 
Taking the sum of equations in (QQ) gives 

which gives a constant net current in the mean field approximation J(x) = Jo within each 
section in the steady state. A positive (negative) net current indicates a net current towards 
the plus (minus) end. 

A. Spatial distribution in the unipolar section 

The system we consider is closed, and so the unipolar sections are half closed, i.e., no 
exit (of plus-type particles) and no injection (of minus-type particles) is possible at the plus 
ends. For the first track, the unipolar section in x 2 < x < 1 is closed at the right end, 
corresponding to the boundary rates a u _ = = 0, which give 

P (l) = l, (7(1) =0. 

We assume the left end x = x<i of the unipolar section is associated with injection/exit rates 
a\ and As discussed in 20J], the unipolar section has a zero net current due to the 
closed boundary at the right end. Thus J = p{x){\ — p(x)) — cr(x)(l — cr(x)) = 0, leading to 
the density relation 

p(x) = 1 — a(x), p(x) = a(x). 
Together with the ODE ([1]), we have two possible solutions - complementary density 

p{x) = l-a(x) = 1-n + Qx (3) 

and equal density p(x) = a(x) = C. The equal density may appear in the bulk of the 
unipolar while the complementary density starts from the plus end and extends toward the 
interior of the unipolar section. 

10 



1 st track 



p! 



© 



p: a; 



t4® 



+ 

p: 



FIG. 2: One of the tracks for ASEP with antipolar bundling described in Figure [T] is mapped onto 
two unbundled bidirectional ASEPs coupled by boundary conditions as shown here. This shows 
the sections of the first track. The boundary rates a" = fft = as the system in Figure [1] is closed. 



As discussed in 
case reads as C = 



20], for low injection and high exit rates, the constant in the equal density 



a 



the density profile of plus-type particles exhibits a shock 



between equal density and complementary density regimes while the density of minus-type 
particles is continuous except near the left end. We refer this as an SL phase; the first letter 
means for the plus-type particles while the second one for the minus-type particles; S stands 
for shock while L stands for low (less than a half). The shock location x s is determined by 
matching lim^^- p(x) = lim a ,_ >a .+ 1 — p{x). Combining the complementary solution ([3]) and 
equal density a" gives x s = 1 — A limit lim Xs ^. X2 a" (x s ) = Qx\ gives a boundary of 
this SL phase in the space When the shock is driven out of the unipolar section, 

the plus-type particles are in high density (over one half) and minus-type particles are in 
low density; we refer as an HL phase; H stands for high density. The linear density profile 
of plus-type particles in the complementary density decreases from 1 at the plus end and 
propagates inside. When a lower bound of 1/2 for plus-type particles in density is reached 
before x is decreased to X2, the density of 1/2 continues when x is further decreased and the 
maximum unidirectional current occurs in the unipolar section. This is referred as an MM 
phase (M stands for maximal unidirectional current, which is 1/4) and only occurs when 
1 — Q + Qx2 < 1/2, i.e., Qx\ > 1/2. By particle-hole symmetry, analogous to the SL phase, 
we have an HS phase where the equal density is over one half and minus-type particles 
experience a shock. In summary, the phase diagrams of density profiles for the unipolar 
section with general left boundary conditions a" and are shown in Figure [3] for both 
Qx\ < 1/2 and Qxi > 1/2. In the coexistence line a" = we have p{x-i) = al = 1 — cr(x2) 



and note that p(l — ^) 



as m 



al in the complementary density, thus a similar discussion 



14J can be applied in the region [x 2 , 1 — <y- 



11 




FIG. 3: Phase diagram of the density profile in the unipolar section with general left boundary 
conditions off and /3" and closed right boundary condition for both cases Qx% < 1/2 (left panel) 
and £lxi > 1/2 (right panel), where x\ = 1 — %2 = 0.2. Numerical examples of density profiles in 
the unipolar section for each phase are illustrated with indicated boundary conditions (ce",/3"). 
The top three density profiles in the middle column are using = 1 while the bottom density 
profile is using Q = 10. In the density profiles, black dots are for plus-type particles while gray 
dots are for the minus-type. 

B. Analysis of spatial distribution along the entire bundle 

In the following subsections, we aim to analyze the density profiles in the entire symmetric 
bundle by considering two ASEPs on the first track with corresponding boundary rates in 
each section. The boundary rates in one section may depend on those in another section. 

In the uncoupled case (gi = q2 = 0), the stochastic process on each track is exactly the 
same as the two-lane ASEP discussed in [l^] assuming equal direction-change rates. Also 
from the discussion of the unipolar section in Section IIII A[ we can see that the density 
profile on the first track consists of three segments in general; an equal-density segment 
occurs in the middle and connects with the boundaries by complementary-density segments 
on its left and right sides which are continuous at x = x\ and 1 respectively. 

For either q\ or q2 positive, the density profiles can exhibit a variety of phases even if the 
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parameters f2 and G are fixed. In the following, we show that a mean-field approximation 
with appropriate boundary conditions for both unipolar and bipolar sections agrees well 
with numerical simulations for a variety of density profiles. The density profiles in the 
unipolar are discussed in Section III I A\ thus we focus on the bipolar section in order to 
understand the density profile along the entire lattice. In the generic case where qiq 2 > 0, 
the net current in the bipolar section on a single track is not necessary zero. For a small net 
current J, one can show that there are "equal" (approximately equal) and "complementary" 
(approximately complementary) density solutions a ~ p, o 1 — p; see Appendix [A] for 
detailed explanations. Particularly, when neither p nor a is close to one half, we approximate 
the constant net current J = p(l — p) — a(l — a) by 

J = sgn(l - 2p)(p - a) and J = sgn(l - 2p)(p + a - 1) (4) 

for "equal" and "complementary" densities respectively. The solutions of density profiles 
for plus-type particle from (lAip can thus be approximated as 

p = —JQx + C e , p = VLx + C c (5) 

for "equal" and "complementary" densities respectively and the corresponding densities of 
minus-type particles can be approximated as 

cr = p + sgn(2p — 1) J, a = 1 — p — sgn(2p — 1) J (6) 

The constants C e t c \ in ([5]) can be seen as functions of the boundary rates a h ± and f!> h ±. 
These boundary rates are associated with the parameters in the model. Note that for a 
minus- type particle crossing the junction x = X2 from the unipolar to the bipolar section, it 
switches to the other track with rate qi (which contributes to the injection of plus-type of 
particles on the other track in the bipolar section) and keeps on the same track with rate 
1 — qi (which contributes to the injection of minus- type particles). Therefore, by the spatial 
symmetry, we approximate the injection rates by 

a b + = o(x2)qi and a b _ = a(x 2 )(l — qi). (7) 

Moreover, a minus-type particle on the second track switches to the first track with rate q 2 
when crossing the junction at x = x 2 (which contributes to the exit of minus-type particles 
in the bipolar section); a plus-type particle steps forward with rate 1 when crossing the 

13 



junction (which contributes to the exit of plus- type particles); both the minus- and plus- 
type particles share the same target site when they move. Thus, by spatial symmetry, we 
have PX + Pi = (l- p{x 2 )) and fi b _ = q 2 (3 b + . This gives 

# = 1 - P(X2) and P b = ( 8 ) 

1 + 92 1 + 92 

C. Phases for low overall densities 

For a low overall density <C 1/2, the unipolar sections have a small number of particles 
in occupancy, thus the unipolar is in an SL phase where densities of both types of particles 
are equally constant away from the plus end; say a which is associated with a shock location 
a;" = 1 — a/Q in the unipolar section from When both types of particles in the bipolar 
section are in low densities, they are dominated by the injection rates a h ± in (J2J) and we 
approximate ct{x2) by o. Thus, in the bipolar section, the density profile of each type is 
governed by the parameters Q, q\ }2 as well as a. If the end-induced switching rate q 2 = 0, 
then queueing appears at minus ends in a similar manner to the queueing at the plus end 
in the unipolar section. In contrast, for a sufficiently large switching rate q 2 , any queuing 
particles at the minus end are expected to move into the unipolar section of another track. 
This leads to low densities for both types of particles in the bipolar section and we refer as 
an LL-SL phase for the entire system; the two letters before the dash are for the bipolar 
section while the other two follow the dash are for the unipolar, e.g., here LL stands for 
the phase in the bipolar while SL for the unipolar section. The explicit expressions of the 
mean-field approximation for the density profiles can be worked out by using "equal" density 
approximations ([5} and (Q for the bipolar section and considering the conservation of G, i.e., 
® = ^ 1 P 2(i-xi) " — • ^ ne detailed calculation of these expressions are given in Appendix [B] and 
Figure H] (a) shows that the mean- field approximation agrees well with numerical simulations 
for this LL-SL phase. 

In this LL-SL phase, a boundary layer may arise near the junctions. When the boundary 
rate f$ h _ is satisfied, a shock of minus-type particles could form near the minus end. We 
refer this case as an LS-SL phase (i.e., in the bipolar section minus-type particles are in 
shock state and plus-type particles are in low density, while in the unipolar section shock 
forms for the plus type and the density of the minus type is low); see Figure @] (b) as an 
example. In the mean-field approximation, a shock for minus-type particles in the bipolar 
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section stabilizes at x = x b s when 



lim cr(x) = 1 — lim a(x). (9) 



,.6+ 



Similar to the LL-SL phase, the density value a which also gives the boundary rate (3_ in (jSJ) 
when approximating p{x2) by a, can in principle be worked out from the association with pa- 
rameters. Hence, we have approximations using "equal" and "complementary" densities (jHJ) 
and ([6]) for right and left sides of the shock in the bipolar section. The detailed calculations 
are not shown here. Figure H] (b) shows that the mean-field approximation agrees well with 
numerical simulations for this phase. 

The transition between LS-SL and LL-SL phases for low overall density will be that the 
shock in the bipolar section is driven to the junction between sections i.e., x\ — > x\. In other 
words, the "equal" density of minus-type particles a(x) in (1B2I) satisfies \im x ^ Xl a(x) = 
1 — a(xx) = f3 b _, which gives 

Recall that a is related to the overall density by (IB5I) . Figure H] (c,d) shows examples 
of this prediction on the border line between LL-SL and LS-SL phases in the space (91,92) 
against numerical simulations where a shock is identified by over one half density of minus- 
type particles at the last but one site to the minus end. Qualitatively, this prediction agrees 
well with the simulations. 



D. Phases for intermediate overall densities 



For intermediate overall densities, it is clear that the system could be in the LS-SL phase 
and thus the HS-HS phases (via particle- hole "symmetry"). However, there is more variety 
of phases in the density profile along the entire system with different parameters. We have 
not attempted to characterize all the possible phases for intermediate density, but in the 
following we show another two cases in details: one with shocks of both-type particles and 
the other with "smooth connection" of one type in the bipolar section. 
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FIG. 4: (a) and (b) show density profiles in LL-SL and LS-SL phases on the first track with x\ = 
l—x 2 = 0.2 and other parameters indicated. For plus-type (minus-type) particles, black (gray) lines 
show the density solution from mean-field approximations and squares (circles) show the averaged 
densities (over T = 60000 s) from Gillespie simulations, (c) and (d) show the transition between 
LL-SL and LS-SL phases for indicated f2 and G. Lines are from the mean- field approximation (|10|) 
while circles are from numerical simulations - a shock is identified if p]vi+i > 0.5. 

1. SS-HL phase 

Note that a small Q could lead the unipolar section to be in an HL phase which could 
give shocks for both types of particles in the bipolar section. This is referred as an SS-HL 
phase; see Figure [5] for an example. In this phase, the density profiles have "equal" density 
in the middle separating two "complementary" density regimes in the bipolar section; the 
four boundary conditions a b ±1 fi b ± given in (J7|) and (|HJ) are all satisfied. Therefore, we have 

p(x x ) =a\< 1/2, p(x 2 ) = 1-Pl, <r(xi) = 1 - fit, *(x 2 ) = a b _ < 1/2 

If p and a in the bipolar section are not close to one half, then from (j3J) we approximate the 
net current by J = (3 b + — a b _ = a b + — [3 h _. We also approximate p(x2), o(x2) in the boundary 
rates by the limit in the complementary density (j3J) as x — > x 2 , i.e., 1 — p{x 2 ) ~ cr{x 2 ) ~ Xifl. 
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FIG. 5: Density profile on the first track with parameters indicated shows an SS-HL phase. The 
markers for simulation and mean-field approximations are as in Figure 0] (a,b) . 

Thus, the net current is approximated by 

J = P+-a b - = a?+-P b -= (j^-^-l + qijxiQ. (11) 

Similar to the discussion of the LL-SL phase, for the SS-HL phase, the relation between the 
shock locations in the bipolar section and the overall density can in principle be worked out 
and would give approximated solutions of density profiles in the bundle by using fl5]) and (J6j) 
again. Figure [H] (a) shows that the approximated solution from the mean-field agrees well 
with numerical simulations. 

Note that the approximation ffTTj) of the net current suggests that the direction of the 
net current is governed by gx,2 and (1 — qi)(l + q-z) = gives a zero net current. Thus we 
are expected to have constant density in a region of the bipolar section if zero net current is 
satisfied as seen in Section [III A[ Particularly, when the overall density = 1/2, there would 
be a "maximum unidirectional current" region in the bulk where p 1 (x) = a l (x) = 1/2. This 
region can be estimated by identifying degenerated shocks (with zero shock height) for both 
types of particles, which can be approximated as + (1 — qi)xi, 1 — — qiX\]. Note also 
that this "maximum unidirectional current" could also appear with an MM phase in the 
unipolar section. 



2. Phases with "smooth connection 



For a non-zero net current J, the "equal" density solution is a monotonic function in 
position seen from ( lAlj) . Thus, simply by increasing the overall density from a low value, 
the maximum value of the density profiles in the "equal" density region, would increase and 
eventually reach one half. Figure M shows density profiles together with a plot of p vs o 
under different 6; particularly, for 6 = 5/16 and 9 = 9/16, the "equal" density regions 
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contain both high and low densities and one of the types (the plus-type in this example) 
exhibits slower change between low and high density than the other type. Moreover, Figure [7] 
compares the density profiles in different system size where the density profile in the middle 
is more shock-like for plus-type particles and remains almost unchanged for the other type. 
These suggest that in the bipolar section p is smoothly increasing with x in the "equal" 
density region while a has a shock profile separating low and high densities in the limiting 
system. We call this phase with smooth connection between low and high density for one 
type of particles as an SC (or a CS) phase (ignoring the existence of shocks for the type 
which has smooth connection in the bipolar section) depending on which type of particles 
exhibits "smooth connection"; the letter 'S' stands for shock while the letter 'C stands for 
connection. Whether the density of plus-type or minus-type particles smoothly connects low 
and high density is related to the sign of the net current J; a positive J is associated with 
plus-type particles having smooth connection and vice versa. 

When in a CS (or an SC) phase in the bipolar section, the phase in the unipolar section 
can be different; see Figure [6] and Figure [8] where the unipolar section can be in an HS, 
HL, SL or SS phase depending on the parameters. Note that the SS phase in the unipolar 
occurs when the boundary rate a" = in contrast to the polynomial function of density 
profiles shown in I4J] with symmetric open boundary condition, both density profiles exhibit 
shocks here due to the overall density being fixed; this is consistent with the discussion in 



14j. Furthermore, the "smooth connection" phases are generic in the parameter space qi^] 



seen from Figure [9] where each frame represents the density profiles for p(x) (left panel) 
and cr(x) (right panel) via color for fixed q± and changing q2 by every 0.1 between and 1. 
The parameter q\ e {0,0.2,0.5,0.8, 1} is increased from top frame to the bottom frame in 
both panels. In addition, there are both SC and CS phases in the parameter space ((71,(72)- 
For instance, plus-type particles exhibit a smooth connection for high qi while minus-type 
particles exhibit for low q\. 

The density profiles in "smooth connection" phases can also be understood by mean- 
field approximations. When using approximation (lAlj) for densities in the bipolar section, 
the constants can in principle be worked out by considering the fact of the overall density 
conservation - similar to the previous discussion on the LL-SL phase. However, due to the 
non-linearity in the expression flAip . we examine the mean-field approximation for the CS 
(or SC) phase by choosing appropriate J and boundary conditions to satisfy the overall 
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0.2 < x < 1 



FIG. 6: Left: density profiles on the first track for four different overall densities as indicated; 
dark for plus-type particles while gray for minus-type particles. Right: plot of p vs a from the 
corresponding left panel; the solid lines show the relation J = p(l — p) — a(l — a) where J is the 
constant net current. Deviation of the dots from the solid lines are probably due to finite size or 
boundary layer effects. Other parameters are O = 5, q% = 0.8, q2 = 0.2. 




0.5 0.6 0.7 

x\ < x < 1 



FIG. 7: Density profiles on the first track for system sizes N = 300 ('+') and N = 1000 ('.'). The 
black is for p while gray is for a in both system sizes. In the middle of the lattice, the shock width 
for p is shortened in the larger system size; whereas a is almost unchanged in two system sizes. 
Other parameters are f2 := u> x N = 5, q\ = 0, qi = 0.5, O = 1/2. 

density G and the densities in the unipolar section. Take Figure [8] (a) as an example where 
minus-type particles smoothly connect low and high densities and also form a shock in the 
bipolar section, and in the unipolar section densities are in an HS phase with an equal high 
density which is constant away from the plus end, say a. In the CS-HS phase, given a net 
current J and a with boundary conditions ot b + = aq 1 and 1 — qi), we would have the 
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FIG. 8: Density profiles on the first track with indicated parameters showing "smooth connection" 
phases for intermediate overall densities. A CS phase in the bipolar section indicates that minus- 
type particles exhibit a smooth connection through density 1/2 at a point where plus- type particles 
exhibit a shock, while for an SC phase the plus-type particles exhibit the smooth connection. The 
markers for simulation and mean-field approximations are as in Figure U] (a,b) . 



"equal" and "complementary" densities by approximation flAlj) in the bipolar section and 
would also have the equal and complementary densities from Section IIII Al in the unipolar 
section. The parameters J and a are chosen in order to match the shock condition for the 
minus-type particle in the bipolar section as in ^ and the overall density G = ^rn~rj ^ — ■ 
Particularly, in this example, the mean-field solution with a = 0.753 and J = —0.065 satisfy 
the shock condition and the overall density; Figure [8] (a) shows the agreements between the 
mean-field approximation with chosen parameters and the simulation. Similar comparisons 
for "smooth connection" phases with an HL, SL or SS phase in the unipolar section are 
shown in Figured] (b-d). 



E. Phases for high overall densities 

When the system is in a high overall density, the vacancies are in low overall density. 
From the discussion of particle-hole "symmetry" in Section [III we can see that an HH-HS 
and SH-HS phase could appear. Moreover, note that such "symmetry" may break when qi 
is large, in which case we find another phase where in both sections plus-type particles are in 
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FIG. 9: A variety of phases arise in the bundle for fixed = 7/16 and $7 = 5 when varying the 
switching rates q\%. Each frame in the left and right panel represents densities of plus- and minus- 
type particles for fixed q\ when varying 52 by every 0.1 between and 1. The fixed parameter 
qi € {0,0.2,0.5,0.8, 1} is increased from top to bottom frame in both panels. 
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FIG. 10: Density profiles on the first track with indicated parameters shows an HS-HS phase. The 
markers for simulation and mean-field approximations are as in Figure 0] (a,b) . 

high density while minus-type particles exhibit a shock, and refere as an HS-HS phase; see 
Figure ITOl In this phase the boundary rates a b _,p± are satisfied and Figure fTUI shows that 
the mean-field approximation agrees with numerical simulations on the density profiles. 



IV. DOMINANCE OF PARTICLES IN THE TRANSPORT 



Both types of particles are essential for bidirectional transport along the bundle of tracks 
in the unipolar section. However, it is possible to have only one type of particles take the 
bidirectional transport across the bipolar section by track switching events and indeed in 
vivo experiments of EEs suggest that kinesin-3 is the main motor for long-range EE motility 
across an antipolar MT bundle 21] . In the following, we discuss the contribution of each 
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type of particles to the transport in terms of occupancy as well as current in our model. 



A. Dominance in occupancy 

For the transport in the entire bundle, one way to quantify the dominance in occupancy 
is to estimate the fraction of plus-type particles in the entire bundle, < F + < 1 where 

_ 6+ = j l Xl p\x)dx 
+ ' 26(1 -xO 

for a symmetric lattice x\ + x 2 = 1- Recall that 6 is an overall density of particles in the 
bundle which is preserved during the transport and is considered as a parameter. It is clear 
that F + = 1/2 in the case qi = q 2 = 0. Note that track switching is accompanied with a 
change of type from minus to plus, thus increasing either qi or q 2 will have the potential 
to increase the number of plus- type particles and thus increase the fraction F + . If F + is 
close to 1 then plus-type particles are in significant dominance in terms of occupancy in the 
entire bundle. For Q = and qi >2 > together with a sufficiently low overall density 0, the 
maximum F + — 1 can be achieved. 

When considering the fraction of plus-type particles within each section, it is clear that in 
the unipolar section there are more plus-type than minus-type particles for any parameters. 
In contrast, the fraction in the bipolar section 

Fb . = C p^ dx 

+ fx? P( x ) + <*(x)dx 
is less easy to estimate for general overall density 0. 

For low overall density where the system is in either LL-SL or LS-SL phase, the mean-field 
approximation predicts the corresponding density profiles well. Thus, for sufficiently large 
q 2 (i-e., larger than the critical value satisfying ffTOj) ) which gives low densities for both types 
of particles in the bipolar section, both fractions F + and i 7 ^ can be expressed using the 
approximated density profiles, ignoring boundary layers. The fraction in the entire bundle 
is 

_ aQ(x 2 - xi) 2 /2 + aq^Xi - xQ + 1 - a 2 /(2Q) + (xi - a/fl)a 
+ ~ a/n + (1 - 2a/n)a 1 ' 

where a is associated with as in (IB5|) . and the fraction in the bipolar section is 

F b ^- + 2qi ~ 1 (13) 
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FIG. 11: The fraction of plus-type particles in the bipolar section. Left: lines are from approx- 
imated equation (|13p for Q € {0.1,1,2,5,10} with bolder lines for larger Q while circles are for 
simulations with larger circle for larger f2. Other parameters are using O = 1/16 and q<i = 1. 
Right panels: to examine the independence of this fraction on the overall density and switching 
rate q 2 if sufficiently large; cross and circle are for = 3/16 and 1/16 respectively in both up and 
down panels. Simulated data presented are using q\ = 0.2 - quantitatively similar figures can be 
obtained for other q\. 

Both fractions show independence of q 2 and linear dependence on q±. Comparatively, the 
latter expression for the fraction in the bipolar section is a simpler function of Q and q± 
for fixed X\p- We can see that this expression agrees well with simulations from Figure fTTl 
Moreover, from this approximation, qi = 1/2 gives approximately equal contribution in the 
bipolar section for each type of particles. This agrees with a zero net current in (IB3[) which 
leads to equal density in an LL phase. Moreover, the sign of q\ — 1/2 determines which 
type of particles are in dominance and by decreasing Q, the dominance can be enhanced. 
Particularly, for a small Q 

lim FUqi) « 1, lim FUqi) « 0. 

qi-tl x gi-s>0 T 

That is to say, for sufficiently low Q, plus- type particles are in significant dominance in the 
bipolar section for sufficiently large switching rate qi, in contrast, minus-type particles are 
in significant dominance for sufficiently small q\. 

For an intermediate overall density, the plus-type and/or minus-type particles form shocks 
near the junctions between sections. In contrast to low overall densities where is highly 
dependent on q\ and weakly dependent on q 2 , in the intermediate overall densities, the 
fraction in the bipolar section F^_ shows relatively weak dependence on both switching rates 
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FIG. 12: The fraction of plus-type particles (left panel) and the ratio of average current for plus- 
type particles R b + (right panel) in the bipolar section as varying q\ (marked as *) or 52 (marked as 
+) for Q = 5 and = 7/16. Simulated data presented are using q\ = 0.8 (#2 = 0.1) when varying 
92(91) - qualitatively similar figures can be obtained for other values of q\ t 2- 

qi t 2- The left panel in Figure H~2l shows as an example with an overall density G = 7/16 of 
how the fraction in occupancy varies changes with the rates qip- 



B. Dominance in current 



The directed currents J± are also important quantities to characterize the transport 
capacity - a large current indicates efficient transport. In contrast to the standard unidirec- 
tional ASEP on a single lane (where both the density and current are constant along the 
bulk of the lane), in our model, neither the density nor the directed current of each type of 
particles on a single track is constant; see Figure [T3l for examples. Thus, we consider the 
average unidirectional current 

fyJ±(x)dx 

\J±)x '■— — : 

1 — X\ 

and define an overall current (J) x '■= (J+)x + (J-)x- The dominance in current of particles 
can be investigated by looking at the ratio of average current for plus-type particles: 

D _ (J+)x _ (J+)x 

+ - (J) ffi (J + ) x + (JJ) X 

We simultaneously define the ratio of average currents within each section. In the unipolar 
sections, plus- and minus-type particles equally contribute to the average currents as the 
net current is zero. Thus the ratio in the entire bundle R + depends only on the ratio within 
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FIG. 13: Numerical current profiles (bottom) with corresponding density profiles (top) for overall 
density © = 3/16 (left) and © = 9/16 (right). Other parameters use Q = 5, qi = 0.8, q2 = 0.4. 
Dark dots are for plus-type particles while gray dots are for the minus-type. 



the bipolar section, 



Sxi J +i X ) dx 



£(J + (x) + J-(x))dx- 
When the density profiles in the bipolar section are in the LL phase (p, a <C l/2) ; the 
dominance in current is equivalent to that in occupancy, ~ F? (seen from Figure [T51 (left 
panel)). However, if shocks form in the bipolar section, higher densities (above one half) can 
give a lower current (seen from Figure fT3l (right panel)). Thus dominance in occupancy does 
not imply dominance in current; see Figure fT2l as an example. Moreover, by comparing the 
two panels in Figure [12], shows a larger range than i 7 ^ on changing parameters. Thus, 
for an intermediate overall density, one could say that current plays a more important role 
than occupancy in determining transport of particles. In addition, we can see that the ratio 
in current R b + shows a larger range on changing q\ than on changing q2 in this example. 
This suggests that q\ is more important than q2 for the dominance. 



V. DISCUSSION 



In this paper, we introduce an ASEP-type model to describe the bidirectional motility of 
particles on an oriented bundle of tracks. This aims to model the motion of cargos/motor 
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complexes undergoing transport along an antipolar MT bundle within a cell. Our model 
is certainly a great simplification of cell transport processes. It is parameterized by the 
particle turning rate Q (the inverse of the run length), the overall density O (the proportion 
of sites that are occupied by particles) and the obstacle/end-induced switching rates (i.e. 
the rate at which particles switch MTs at the junctions between the unipolar and bipolar 
sections). We use numerical simulations and mean-field approximation to investigate the 
dependence of the stationary density profiles within the bundle on these parameters. 

We observe that, as expected, the switching rates gi >2 have a major effect on the distribu- 
tion of particles along the bundle. Although we have not fully explored the dependence of 
phases on parameters we highlight below a number of interesting features about the system. 
Even for low overall density G, particles (cargo/motor complexes) can queue to form accu- 
mulations at minus ends in addition to any plus-end queuing observed in simpler situations 
[20I . The critical value of the end- induced switching rate g 2 above which shocks are formed 
at minus ends is investigated in Section fill CI For intermediate overall density G, we find a 
variety of phases, including a new "smooth connection" phase where in the limit, the density 
profile of one type of particle on a track smoothly passes through one half while the other 
type displays a shock. We discuss a variety of other phases and have investigated the role of 
the switching rates gi )2 in determining which motor is dominant in the transport. For low 
overall densities, the obstacle-induced switching rate q\ influence the fraction of particles of 
one type within the bipolar section in an approximately linear way. For high overall density 
G, our study shows that although the fraction of occupancy by different type particles does 
not vary much with rates gi j2 , the fraction that actually contributes to the current (i.e., 
transport) from different types may vary much more. 

Our model is inspired by in vivo experimental observations, although the model has 
been simplified in many ways as we now discuss. For example, it is probably that hopping 
between MTs is not restricted to the ends of the bipolar section, but may occur throughout 
the bipolar section [2JJ. This would allow more possibilities for transition events, and the 
possibility of plus-type particles (in addition to minus-type particles) switching MTs. A 
switch in transport direction can result from hopping between MTs and the activity of 
kinesin-3 alone or it could be a consequence of dynein binding to the cargo, which in case of 
EEs was shown to override kinesin-3 activity [19|. Our model also assumes only two lanes 
in each motility track. Moreover, it is known that numerous MTs form a bundle 21|, each 
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consisting of 13 protofilaments [17] . Thus, many more tracks might support bidirectional 
motility of the cargo/motor complexes and allowing opposite-directed particles moving on 



the same protofilament wil 
plus ends, as discussed in 



certainly give new effects such as increasing the jamming at MT 



161 ] . The complex geometry of the bundle will also contribute to 
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new effects of cooperative transport, already considered for unidirectional transport in 

Notwithstanding these simplifications, we suggest the model could be useful in a number 
of ways, especially when improved measurement of transition rates is possible in vivo, and 
we highlight some of these below: 

1. For low overall density of particles, it is possible to have an accumulation of particles 
at minus ends of MTs. As in vivo experiments so far show no obvious accumulation 
of early endosome cargos throughout the entire cell {21], if this model is accurate then 
it suggests the end-induced switching rate must be high in vivo. There are clearly 
other possible explanations due to features not included in the model such as a high 
turning rate of minus-directed organelles at minus ends (which might be due to an 
accumulation of motors at the minus end taking dynein to plus ends) could avoid 
such an accumulation of EEs; this would be a similar mechanism to that suggested 
in 
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201 ] where an accumulation of dynein motors increases the turning rate of EEs 



near plus-ends and so avoids accumulations of EEs. 

2. For low overall density of cargos, the fraction of occupancy can be used to show how 
the proportion of EEs carried by kinesin-3 among all EEs (carried by either kinesin-3 
or dynein) varies with the parameters in Section HVl These estimates could in principle 
be explored experimentally to test the modelling assumptions used to describe in vivo 
bidirectional transport. 

Further experimental work will improve our qualitative and quantitative knowledge of the 
details of transport process in the living cell, though the consequences of individual cargo 
behaviour on coordinated transport are not always clear. By using the model discussed in 
this paper, and developments thereof, there is a platform with which one can address various 
open questions about the nature and function of the coordinated systems involved in cell 
transport 
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Appendix A: General density profiles in the bipolar section 



In the bipolar section, the net current J (as discussed in Section IHip given by J = 
p(l — p) — er(l — a) is constant when the system is in a statistically stationary state. This 
density-current relation gives 



2 2 \\(2p-iy ) 2 

which implies the "equal" and "complementary" density solutions: 



Thus, the first order ODE (PQ) from the mean field approximation reads 



„ dp Q I 4J 

= -r- + — ±17 — -f- 1 - 1 

dx 2 y Y (2/o - l) 2 
with a general solution of "equal" or "complementary" density 

4J N 3/2 



I I {2 p-iy ) ± 1 j( 2 P- 1 ) 3 = Tl2Jttx + C (Al) 
for any non-zero net current J on taking options in "±" and "=f" . 

Appendix B: Mean-field analysis for an LL-SL phase 

When both types of particles in the bipolar section are in low densities, they are domi- 
nated by the injection rates a b ± which are approximated by 

a b + = aqi and a b _ = a{\ — q{). (Bl) 

where a is the equal density in the unipolar section. Hence from (J^J) and dHJ), the densities 
on the first track in the bipolar section are approximated by 

p = -JQ(x-x 1 ) + a b + , a = -m(x-x 2 ) + a b _. (B2) 

Together with the approximation J = p — a, we find 

= a b + _ a b_ - a{2qi _ 1} 

Q(x 2 -x 1 ) + l Q(x 2 -x 1 ) + V 1 ' 
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which gives a zero net current if qi = 1/2. The density expressions (1B2[) together with the 
symmetry p 1 {x) = p 2 (l — x) and o l (x) = <r 2 (l — x) give the following 

p 1 (x) + a 2 (x) = p 2 (x) + o 1 (x) — a, x G x 2 ) ■ 

Thus 

£ P( x ) + °{x)dx _ g/n + (1 -2a/n)a 

2{l-x l ) ~ 2(1 -xO ' ( ] 

That is, given a low overall density 9, a sufficiently large q 2 and assuming densities are low 
in the bipolar section, we have 



Substituting f lB5j) . (JB3J) and (1B1I) into (1B2j) gives the analytical approximation of the density 
profiles in the bipolar section in an LL-SL phase. This, together with density profiles in the 
unipolar sections, gives the density profiles in the entire bundle. 
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